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We find that coherent multiple scatterings of the partonic fluctuations of the virtual 
photon with the nucleons inside the nucleus result in the large nuclear shadowing effect 
for parton distributions in nuclei. We predict the gluon shadowing effect at small x and 
for Q 2 > 3GeV 2 . 
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The theoretical studies of semihard processes play an increasingly important role in 
ultrarelativistic heavy-ion collisions at > 200 GeV in describing the global collision 
features such as particle multiplicities and transverse energy distributions. These semihard 
processes can be reliably calculated in the framework of perturbative QCD and they are 
crucial for determining the initial condition for the possible formation of a quark-gluon 
plasma. Assuming the validity of the factorization theorem in perturbation theory, it is 
essential to know the parton distributions in nuclei in order to compute these processes. 

The attenuation of quark density in a nucleus has been firmly established experimen- 
tally at CERN M and Fermilab in the region of small x in deeply-inelastic lepton 
scatterings (DIS) off nuclei. The data, taken over a wide kinematic range, 10 -5 < x < 0.1 
and 0.05 GeV 2 < Q 2 < 100 GeV 2 , show a systematic reduction of nuclear structure func- 
tion F^(x, Q 2 )/A with respect to the free nucleon structure function F^(x,Q 2 ). There 
are also some indications of nuclear gluon shadowing from the analysis of J/\l/ suppres- 
sion in hadron-nucleus experiments |J but the extraction of nuclear gluon density is not 
unambiguous since it involves the evaluation of the initial parton energy loss and final 
state interactions ||. 

At low Q 2 , in DIS, the interaction of the virtual photon with the nucleons in the rest 
frame of the target is most naturally described by a vector-meson-dominance (VMD) 
model [[J. At Q 2 > 1 ~ 2 GeV 2 , the virtual photon can penetrate the nucleon and probe 
the partonic degrees of freedom where a partonic interpretation based on perturbative 
QCD is most relevant in the infinite momentum frame. In the target rest frame, the virtual 
photon interacts with nucleons via its quark- ant iquark pair (qq) color-singlet fluctuation 
||. If the coherence length of the virtual photon is larger than the distance between 
nucleons in a nucleus, l c > Rnn, the qq configuration interacts coherently with fraction 
of the nucleons, while for l c > Ra (i.e. x < 10 -2 ), it interacts coherently with all the 
nucleons, with a cross section given by the color transparency mechanism for a point-like 
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color-singlet configuration, a qq N = ^-r 2 a s x' gr>LA(x' ', l/ r t ) H- This results in a reduction 
of the total cross section and consequently attenuation of the parton distributions in 
nuclei. 

In the Glauber-Gribov multiple-scattering theory |7J nuclear collision is a succession of 
collisions of the probe with individual nucleons within nucleus. A partonic system (h), 
being a qq or a gg fluctuation, can scatter coherently from several or all nucleons during 
its passage through the target nucleus. The interference between the multiple scattering 
amplitudes causes a reduction of the hA cross section compared to the naive scaling result 
of A times the respective hN cross section, the origin of the nuclear shadowing. The total 
hA cross section is given by 



J d 2 b2 [l - e -«™^(b)/2] = 2nR 2 A [7 + lnKh + Ei{Kh)] } (1) 

where Kh = Auhn / {^Ra) ls an impact parameter averaged effective number of scat- 
terings. For small value of Kh, &hA 27rR A Kh — A<ThN, the total hA cross section is 
proportional to A. In the limit Kh — > oo, the destructive interference between multiple 
scattering amplitudes reduces the cross section, OhA — > ^Ra(i + ^ nK h)- Namely, the ef- 
fective number of scatterings is large and the total cross section approaches the geometric 
limit 2nR A , a surface term which varies roughly as A 2 / 3 . 
In the Glauber-Gribov eikonal approximation, 

a(>y*A) = / d 2 b f 1 dz [ rf 2 r|^(^,r)| 2 2(l - e -^n{v,z)T A {b)/2^ ^ 



where \ip(z, r)| is the photon wave function |T0| and z is the fraction of the energy carried 



by the quark (antiquark). The nuclear cross section is therefore reduced when compared 
to the simple addition of free nucleon cross sections. 

At small x, the structure function of a nucleus, F^(x, Q 2 ) can be obtained from <j{^f* A). 
Substituting integration over (z, r t ) to (x',Q' 2 ) in (2), one obtains for (sea)quarks 

xf A (x,Q 2 ) = xf A (x,Ql) + -^x — dQ' 2 + ln(K q ) + E 1 (K q )} , (3) 

o7T Jx X J Qq 

where K q (x,Q 2 ) = 3 ^Q 2 <y s {Q 2 )xg^ hA {x, Q 2 ). We find our result for F^(x,Q 2 ) to be in 

excellent agreement with the recent HERA data for a broad range of x and Q 2 

We parametrize the initial shadowing ratio Rq{x) at Qq = 0.4 GeV 2 using experimental 
data and we calculate the nuclear structure function at the measured (Q 2 ) values at 
different x values. The results for 40 Ca and 208 Pb are shown in Fig. 1. We note that at 
large Q 2 the nuclear shadowing effect is reduced but does not diminish, i.e. it is not a 
higher twist effect. This can be understood as the interplay between the perturbative and 
the non-perturbative shadowing mechanisms, which is also evident in the x-dependence 
of R q . The apparent flatness of the shadowing ratio at low Q 2 in the small-x region 
is altered by the perturbative evolution. This is due to the singular behavior of xg^ LA 
as x —> at large Q 2 leading to the strong x-dependence of the effective number of 
scatterings, K q (x). 

Similarly, gluon distribution in a nucleus is given by 

2R 2 A r 1 dx' r® 2 

xg A (x, Q 2 ) = xg A (x, Q 2 ) + \ — dQ' 2 [ 7 + ln(« fl r) + E 1 {k 9 t)\ , (4) 

7T Jx X JQq 
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Fig. 1 

The essential difference between the quark and gluon cases is different splitting functions 
and the effective number of scatterings, n g = 9n g /4: due to different color representations 
that the quark and the gluon belong to. These two effects combined result in the 12 
times faster increase of the gluon density with Q 2 than in the case of the sea quarks in 
the region of small x. The two important effects which make the gluon shadowing quite 
different from the quark shadowing are the stronger scaling violation in the semihard scale 
region and a larger perturbative shadowing effect. This can be seen by considering the 
shadowing ratio 



R g (x,Q 2 



xg A (x, Q 2 ) xg N (x, Q 2 )R°Jx) + Axg A (x; Q 2 , Q 



Axg N (x,Q 2 ) 



xg N (x, Ql) + Axg N (x; Q 2 , Q\ 



(5) 



where R g {x) is the initial shadowing ratio at Ql and Axg(x; Q 2 , Ql) is the change of the 
gluon distribution as the scale changes from Ql to Q 2 . The strong scaling violation due to 
a larger n g at small x causes AxgN^x; Q 2 , Ql) ^> xgN^x, Ql) as Q 2 is greater than 1 ~ 2 
GeV 2 . As seen in Fig. 2(a) for Q 2 > 3GeV 2 the dependence of R g (x,Q 2 ) on the initial 
condition R° g {x) diminishes and the perturbative shadowing mechanism takes over. 

The ^-dependence of the gluon shadowing can also be predicted as long as Q 2 > 3GeV 2 
where the influence of the initial condition is minimal. The shape of the distribution is 
quite robust in the small-x region regardless of what initial conditions one may choose. 
Due to the perturbative nature of the shadowing, these distributions do not exhibit a 
saturation as x decreases. Furthermore since the shadowing is a non-linear effect in the 
effective number of scatterings, the impact parameter dependent shadowing ratio cannot 
be factorized into a product of an average shadowing ratio and the nuclear thickness 
function. Our results are presented in Fig. 2(b). 

In summary, the nuclear shadowing phenomenon is a consequence of the parton co- 
herent multiple scatterings. While the quark density shadowing arises from an interplay 
between "soft" physics and the semihard QCD process, the gluon shadowing is largely 
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driven by a perturbative shadowing mechanism due to the strong scaling violation in the 
small-x region. The gluon shadowing is thus a robust phenomenon at large Q 2 and can 
be unambiguously predicted by perturbative QCD. The strong scaling violation of the 
nucleon structure function in the semihard momentum transfer region at small x can be 
reliably described by perturbative QCD and is a central key to the understanding of the 
scale dependence of the nuclear shadowing effect. The impact parameter dependence of 
gluon shadowing is a non-linear effect in the nuclear thickness function. It is important to 
correctly incorporate the impact parameter dependence of the nuclear structure function 
when one calculates the QCD processes in the central nuclear collisions. 
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